Understanding how transmembrane domains regulate interactions between human BST-2 and the SARS-CoV-2 accessory protein ORF7a.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID, replicates at intracellular membranes. Bone marrow stromal antigen 2 (BST-2; tetherin) is an antiviral response protein that inhibits transport of viral particles after budding within infected cells. RNA viruses such as SARS-CoV-2 use various strategies to disable BST-2, including use of transmembrane 'accessory' proteins that interfere with BST-2 oligomerization. ORF7a is a small, transmembrane protein present in SARS-CoV-2 shown previously to alter BST-2 glycosylation and function. In this study, we investigated the structural basis for BST-2 ORF7a interactions, with a particular focus on transmembrane and juxtamembrane interactions. Our results indicate that transmembrane domains play an important role in BST-2 ORF7a interactions and mutations to the transmembrane domain of BST-2 can alter these interactions, particularly single-nucleotide polymorphisms in BST-2 that result in mutations such as I28S. Using molecular dynamics simulations, we identified specific interfaces and interactions between BST-2 and ORF7a to develop a structural basis for the transmembrane interactions. Differences in glycosylation are observed for BST-2 transmembrane mutants interacting with ORF7a, consistent with the idea that transmembrane domains play a key role in their heterooligomerization. Overall, our results indicate that ORF7a transmembrane domain interactions play a key role along with extracellular and juxtamembrane domains in modulating BST-2 function.

Predicting the Assembly of the Transmembrane Domains of Viral Channel Forming Proteins and Peptide Drug Screening Using a Docking Approach.

A de novo assembly algorithm is provided to propose the assembly of bitopic transmembrane domains (TMDs) of membrane proteins. The algorithm is probed using, in particular, viral channel forming proteins (VCPs) such as M2 of influenza A virus, E protein of severe acute respiratory syndrome corona virus (SARS-CoV), 6K of Chikungunya virus (CHIKV), SH of human respiratory syncytial virus (hRSV), and Vpu of human immunodeficiency virus type 2 (HIV-2). The generation of the structures is based on screening a 7-dimensional space. Assembly of the TMDs can be achieved either by simultaneously docking the individual TMDs or via a sequential docking. Scoring based on estimated binding energies (EBEs) of the oligomeric structures is obtained by the tilt to decipher the handedness of the bundles. The bundles match especially well for all-atom models of M2 referring to an experimentally reported tetrameric bundle. Docking of helical poly-peptides to experimental structures of M2 and E protein identifies improving EBEs for positively charged (K,R,H) and aromatic amino acids (F,Y,W). Data are improved when using polypeptides for which the coordinates of the amino acids are adapted to the Cα coordinates of the respective experimentally derived structures of the TMDs of the target proteins.

A Uniquely Stable Trimeric Model of SARS-CoV-2 Spike Transmembrane Domain.

Understanding fusion mechanisms employed by SARS-CoV-2 spike protein entails realistic transmembrane domain (TMD) models, while no reliable approaches towards predicting the 3D structure of transmembrane (TM) trimers exist. Here, we propose a comprehensive computational framework to model the spike TMD only based on its primary structure. We performed amino acid sequence pattern matching and compared the molecular hydrophobicity potential (MHP) distribution on the helix surface against TM homotrimers with known 3D structures and selected an appropriate template for homology modeling. We then iteratively built a model of spike TMD, adjusting "dynamic MHP portraits" and residue variability motifs. The stability of this model, with and without palmitoyl modifications downstream of the TMD, and several alternative configurations (including a recent NMR structure), was tested in all-atom molecular dynamics simulations in a POPC bilayer mimicking the viral envelope. Our model demonstrated unique stability under the conditions applied and conforms to known basic principles of TM helix packing. The original computational framework looks promising and could potentially be employed in the construction of 3D models of TM trimers for a wide range of membrane proteins.

Ni2+ Catalyzed Cleavage of TrpLE-Fused Small Transmembrane Peptides.

In addition to a membrane anchor, the transmembrane domain (TMD) of single-pass transmembrane proteins (SPTMPs) recently has shown essential roles in the cross-membrane activity or receptor assembly/clustering. However, these small TMD peptides are generally hydrophobic and dynamic, difficult to be expressed and purified. Here, we have integrated the power of TrpLE fusion protein and a sequence-specific nickel-assisted cleavage (SNAC)-tag to produce small TMD peptides in a highly efficient way under mild conditions, which uses Ni2+ as the cleavage reagent, avoiding the usage of toxic cyanogen bromide (CNBr). Furthermore, this method simplifies the downstream protein purification and reconstitution. Two representative TMDs, including the Spike-TMD from severe acute respiratory syndrome coronavirus 2 (SARS2), were successfully produced with high-quality nuclear magnetic resonance (NMR) spectra. Therefore, our study provides a more efficient and practical approach for general structural characterization of the small TM proteins.

Secondary Structures of the Transmembrane Domain of SARS-CoV-2 Spike Protein in Detergent Micelles.

Spike protein of SARS-CoV-2 contains a single-span transmembrane (TM) domain and plays roles in receptor binding, viral attachment and viral entry to the host cells. The TM domain of spike protein is critical for viral infectivity. Herein, the TM domain of spike protein of SARS-CoV-2 was reconstituted in detergent micelles and subjected to structural analysis using solution NMR spectroscopy. The results demonstrate that the TM domain of the protein forms a helical structure in detergent micelles. An unstructured linker is identified between the TM helix and heptapeptide repeat 2 region. The linker is due to the proline residue at position 1213. Side chains of the three tryptophan residues preceding to and within the TM helix important for the function of S-protein might adopt multiple conformations which may be critical for their function. The side chain of W1212 was shown to be exposed to solvent and the side chains of residues W1214 and W1217 are buried in micelles. Relaxation study shows that the TM helix is rigid in solution while several residues have exchanges. The secondary structure and dynamics of the TM domain in this study provide insights into the function of the TM domain of spike protein.

Analysis of SARS-COV2 spike protein variants among Iraqi isolates.

The ongoing pandemic of COVID-19 caused by the SARS-COV2 virus has triggered millions of deaths around the globe. Emerging several variants of the virus with increased transmissibility, the severity of disease, and the ability of the virus to escape from the immune system has a cause for concerns. Here, we compared the spike protein sequence of 91 human SARS CoV2 strains of Iraq to the first reported sequence of SARS-CoV2 isolate from Wuhan Hu-1/China. The strains were isolated between June 2020 and March 2021. Twenty-two distinct mutations were identified within the spike protein regions which were: L5F, L18F, T19R, S151T, G181A, A222V, A348S, L452 (Q or M), T478K, N501Y, A520S, A522V, A570D, S605A, D614G, Q675H, N679K, P681H, T716I, S982A, A1020S, D1118H. The most frequently mutations occurred at the D614G (87/91), followed by S982A (50/91), and A570D (48/91), respectively. In addition, a distinct shift was observed in the type of SARS-COV2 variants present in 2020 compared to 2021 isolates. In 2020, B.1.428.1 lineage was appeared to be a dominant variant (85%). However, the diversity of the variants increased in 2021, and the majority (73%) of the isolated were appeared to belong to B.1.1.7 lineage (VOC/alpha variants). To our knowledge, this is the first major genome analysis of SARS-CoV2 in Iraq. The data from this research could provide insights into SARS-CoV2 evolution, and can be potentially used to recognize the effective vaccine against the disease.

Characterizing genomic variants and mutations in SARS-CoV-2 proteins from Indian isolates.

SARS-CoV-2 is mutating and creating divergent variants by altering the composition of essential constituent proteins. Pharmacologically, it is crucial to understand the diverse mechanism of mutations for stable vaccine or anti-viral drug design. Our current study concentrates on all the constituent proteins of 469 SARS-CoV-2 genome samples, derived from Indian patients. However, the study may easily be extended to the samples across the globe. We perform clustering analysis towards identifying unique variants in each of the SARS-CoV-2 proteins. A total of 536 mutated positions within the coding regions of SARS-CoV-2 proteins are detected among the identified variants from Indian isolates. We quantify mutations by focusing on the unique variants of each SARS-CoV-2 protein. We report the average number of mutation per variant, percentage of mutated positions, synonymous and non-synonymous mutations, mutations occurring in three codon positions and so on. Our study reveals the most susceptible six (06) proteins, which are ORF1ab, Spike (S), Nucleocapsid (N), ORF3a, ORF7a, and ORF8. Several non-synonymous substitutions are observed to be unique in different SARS-CoV-2 proteins. A total of 57 possible deleterious amino acid substitutions are predicted, which may impact on the protein functions. Several mutations show a large decrease in protein stability and are observed in putative functional domains of the proteins that might have some role in disease pathogenesis. We observe a good number of physicochemical property change during above deleterious substitutions.

Intracellular autoactivation of TMPRSS11A, an airway epithelial transmembrane serine protease.

Type II transmembrane serine proteases (TTSPs) are a group of enzymes participating in diverse biological processes. Some members of the TTSP family are implicated in viral infection. TMPRSS11A is a TTSP expressed on the surface of airway epithelial cells, which has been shown to cleave and activate spike proteins of the severe acute respiratory syndrome (SARS) and the Middle East respiratory syndrome coronaviruses (CoVs). In this study, we examined the mechanism underlying the activation cleavage of TMPRSS11A that converts the one-chain zymogen to a two-chain enzyme. By expression in human embryonic kidney 293, esophageal EC9706, and lung epithelial A549 and 16HBE cells, Western blotting, and site-directed mutagenesis, we found that the activation cleavage of human TMPRSS11A was mediated by autocatalysis. Moreover, we found that TMPRSS11A activation cleavage occurred before the protein reached the cell surface, as indicated by studies with trypsin digestion to remove cell surface proteins, treatment with cell organelle-disturbing agents to block intracellular protein trafficking, and analysis of a soluble form of TMPRSS11A without the transmembrane domain. We also showed that TMPRSS11A was able to cleave the SARS-CoV-2 spike protein. These results reveal an intracellular autocleavage mechanism in TMPRSS11A zymogen activation, which differs from the extracellular zymogen activation reported in other TTSPs. These findings provide new insights into the diverse mechanisms in regulating TTSP activation.

Mutational insights into the envelope protein of SARS-CoV-2.

The ongoing mutations in the structural proteins of SARS-CoV-2 are the major impediment for prevention and control of the COVID-19 disease. Presently we focused on evolution of the envelope (E) protein, one of the most enigmatic and less studied protein among the four structural proteins (S, E, M and N) associated with multitude of immunopathological functions of SARS-CoV-2. In the present study, we comprehensively analyzed 81,818 high quality E protein sequences of SARS-CoV-2 globally available in the GISAID database as of 20 August 2020. Compared to Wuhan reference strain, our mutational analysis explored only 1.2 % (982/81818) mutant strains undergoing a total of 115 unique amino acid (aa) substitutions in the E protein, highlighting the fact that most (98.8 %) of the E protein of SARS-CoV-2 strains are highly conserved. Moreover, we found 58.77 % (134 of 228) nucleotides (nt) positions of SARS-CoV-2 E gene encountering a total of 176 unique nt-level mutations globally, which may affect the efficacy of real time RT-PCR-based molecular detection of COVID-19. Importantly, higher aa variations observed in the C-terminal domain (CTD) of the E protein, particularly at Ser55-Phe56, Arg69 and the C-terminal end (DLLV: 72-75) may alter the binding of SARS-CoV-2 Envelope protein to tight junction-associated PALS1 and thus could play a key role in COVID-19 pathogenesis. Furthermore, this study revealed the V25A mutation in the transmembrane domain which is a key factor for the homopentameric conformation of E protein. Our analysis also observed a triple cysteine motif harboring mutation (L39M, A41S, A41V, C43F, C43R, C43S, C44Y, N45R) which may hinder the binding of E protein with spike glycoprotein. These results therefore suggest the continuous monitoring of the structural proteins including the envelope protein of SARS-CoV-2 since the number of genome sequences from across the world are continuously increasing.

Modeling of the adsorption of a protein-fragment on kaolinite with potential antiviral activity.

This work aimed at studying the potentiality of interactions between kaolinite surfaces and a protein-fragment (350-370 amino acid units) extracted from the glycoprotein E1 in the transmembrane domain (TMD) of hepatitis C virus capsid. A computational work was performed for locating the potential electrostatic interaction sites between kaolinite aluminol and siloxane surfaces and the residues of this protein-fragment ligand, monitoring the possible conformational changes. This hydrated neutralized kaolinite/protein-fragment system was simulated by means of molecular modeling based on atomistic force fields based on empirical interatomic potentials and molecular dynamic (MD) simulations. The MD calculations indicated that the studied protein-fragment interacted with the kaolinite surfaces with an exothermic process and structural distortions were observed, particularly with the hydrophilic aluminol surface by favorable adsorption energy. The viral units isolation or trapping by the adsorption on the kaolinite nanoparticles producing structural distortion of the peptide ligands could lead to the blockage of the entry on the receptor and hence a lack of viral activity would be produced. Therefore, these findings with the proposed insights could be an useful information for the next experimental and development studies in the area of discovering inhibitors of the global challenged hepatitis and other pathogenic viruses based on the phyllosilicate surface activity. These MD studies can be extended to other viruses like the COVID-19 interacting with silicate minerals surfaces.

Implications for SARS-CoV-2 Vaccine Design: Fusion of Spike Glycoprotein Transmembrane Domain to Receptor-Binding Domain Induces Trimerization.

The ongoing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic presents an urgent need for an effective vaccine. Molecular characterization of SARS-CoV-2 is critical to the development of effective vaccine and therapeutic strategies. In the present study, we show that the fusion of the SARS-CoV-2 spike protein receptor-binding domain to its transmembrane domain is sufficient to mediate trimerization. Our findings may have implications for vaccine development and therapeutic drug design strategies targeting spike trimerization. As global efforts for developing SARS-CoV-2 vaccines are rapidly underway, we believe this observation is an important consideration for identifying crucial epitopes of SARS-CoV-2.